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Abstract 

o 

O 

£NJ ' Using the perturbative QCD amplitudes for B — > tttt and B —* Kir, we have 

Qh| performed an extensive study of the parameter space where the theoretical 

, predictions for the branching ratios are consistent with recent experimental 

data. From this allowed range of parameter space, we predict the mixing 
induced CP asymmetry for B — > ir + ir~ with about 11% uncertainty and 
the other CP asymmetries for B — > tttt, Kit with 40% ~ 47% uncertainty. 
These errors are expected to be reduced as we restrict the parameter space by 
ON \ studying other decay modes and by further improvements in the experimental 



data. 
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I. INTRODUCTION 



o 

o 

23 
sx: 

> 

X 

the weak phase <pi and there is no uncertainty from hadronic matrix elements |T|. B factory 
is expected to yield information, not only on B — > J /ipKg asymmetry but also on various 
other B meson decays. It has been predicted that B — > Kit, tttt decay modes may have 
large CP asymmetries |2],|3[]. While branching ratios for these modes are very sensitive to 
the input parameters, CP asymmetries are expected to be less sensitive. We report here the 
sensitivity of predicted CP asymmetries in the parameter region restricted by experimental 
values for the branching ratios. 



The mixing induced CP asymmetry for B — > J/i/jKs has been shown to depend on only 
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The time dependent CP asymmetry for B°(B ) — ► / transition, where / is CP eigenstate, 
is given by j|: 

(|A/| 2 - 1) cos(AMt) + 2 lm\ f sin(AMt) 

a/(t) " p^TF+i ' {) 

where A/ = (q/p)p(f), and p(/) is defined by the ratio of the decay amplitudes, p(f) = 
A(B — > f)/A(B — > f). For the sake of convenience, we denote the direct CP asymmetry 
proportional to cos(AMt) by a^ r (/) and the mixing induced CP asymmetry proportional 
to sin(AMt) by a mix (f), 

I'M 2 ~~ 1 / p\ 2ImA/ 

adirU) = |^ | 2 - j-, a m i X {J) = | - -. (2) 

In general, a decay amplitude has two kinds of contributions: so-called tree amplitude; and 
penguin amplitude. The decay amplitude for B — > / can be written as 

A(P^/)=£ r T-£ P P, (3) 

where £t,p are Kobayashi-Maskawa(KM) matrix elements, and T, P are decay amplitudes 
with strong final state interaction phases for tree and penguin contributions, respectively 
For example, £r = V^V^ £p = V; b V td for B° -> tt+tt-, and £t = V^V US , £ P = V£V U for 
P° — > K + n~. Defining CP transformation by CP\B) = \B) and CP\f) = T)f\f), the decay 
amplitude for the charge conjugated mode, B — > f, can be written as 

A(B^f)= Vf (C T T-C P P)- (4) 
Thus the general expression for A/ is given as 

where = |£p/£r|> r amp = P/T 1 , and = arg(V^ ) . I n the above expression, 
and arg(£p/£p) have to be invariant under any rotation of quark phases, qj — > e %ai qj. CP 
asymmetry is classified into the following four types depending on the relationship between 
tree and penguin contributions. 

case 1. If £p has the same weak phase as £p, arg(£p/£ T ) = 0, the hadronic matrix elements 
are canceled in A/. Then the direct CP asymmetry vanishes and the mixing induced 
CP asymmetry is strictly given by the weak phase 0, 

adir(f) = 0, a mix (f) = 77/sin(20). (6) 

In the B — > J/tpKs decay mode, this is the case The experimental data 

0-mixi.J '/if>Ks) makes it possible to determine the weak phase 0i = n — &rg(V t * b VtdV c bV* d ) 
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case 2. If the tree contribution is much larger than that of the penguin, \£tT\ 3> |£p-P|, A/ 
is expressed only by KM matrix elements just like in the case 1., A/ = f]fe l2 ^. This 
leads to the same CP asymmetries as shown in eq.(^). When there is no interference 
between tree and penguin contribution, the direct CP asymmetry does vanish, and the 
mixing induced CP asymmetry is directly related to the weak phase 0. 

case 3. If the penguin contribution is much larger than that of the tree, \£pP\ 3> |£t^1, A/ is 
expressed only by the angle between V t * b V t d and £p, 

\ f = Vf e i2ar ^ v ^. (7) 

In this case, the direct CP asymmetry vanishes, and the mixing induced CP asymmetry 
is directly related to some weak phase likewise, 

cidir(f) = 0, a mix (f) = r?/sin(20 / ), (8) 

where 0' = arg(^Vtd^p). In the B — > (pKs, Ksn° decay modes, this is the case. 
(1) B — > 4>Ks decay mode: 

There is no tree contribution in this decay mode, T = 0. The weak phase 0' for this 
decay mode is defined by uf^ in the ds triangle [0], 

0' = argflfclO&Vc.) = (9) 
The CP asymmetry is strictly given by the weak phase uf^, 

a mix (0K s ) = - S m(2uf:). (10) 

This makes it possible to extract the weak phase ujf^ from the measurement of the 
mixing induced CP asymmetry a mix (<j)K ■ Here a>f c s is related to 0i as follows, 



where u s J = arg(V cs V^/V ts V t l) = vr + C(A 2 ). a mix ((pKs) ^ sin2(0i - uj s ^) implies the 
presence of new physics f||§. 

(2) B — > Kg'K decay mode: 

While we realize that it is difficult to measure the time dependence of a Ks7T o at current 
B factories, we will show that B — > Ks^ Q decay mode satisfies |£tT|/|£pP| ~ 2 x 10~ 2 . 
Here £ T = V cd V c * s V* b V us , £ P = V cd V: s V; b V ts , where V cd V* s factor comes from the K°-K° 
mixing. Then, the mixing induced CP asymmetry for B —>■ Ks7T° decay mode is given 



by the weak phase ufj?, 



tc 

ds 



(i^vr u ) = -sin(2<?). (12) 
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Eq.(|i~2|) has been obtained by Ref. ||. But they have assumed that B° — > 7r + 7r~ 
and S° — > K + tt~ decay modes are dominated by the tree and penguin contributions, 
respectively. These assumptions are inconsistent with experiment. Belle Collaboration 
has found the ratio of the branching ratios JLD) , 



BR(B° -> K + 



71 



_ n fin +0.25+0.11 /iq\ 

2 BR(5° - #%°) " U - bU -o-29-o-i6- W 
Neglecting the tree contribution, isospin analysis leads to a conclusion that the above 
ratio is equal to 1. The eq.(|T3|) implies that the tree contribution can not be neglected 
in the B° — > K + tt~ decay mode. 

case 4. If the penguin contribution is comparable with that of the tree, the direct CP asymme- 
try does not necessarily vanish, and the mixing induced CP asymmetry has impurities 
from the penguin or tree contribution. B — > tttt and B° — > K + tt~ fall into this case. 

In B — > Tin decay mode, both KM factors, £t and £p, have the same order of mag- 
nitude, and the tree contribution can interfere with the penguin amplitude. The CP 
asymmetry is parameterized as 

2rsin5sin0 2 M A , 

ddir = 2 ■= — , 14 

1 + r z + 2r cos o cos <p2 



7171 



a dir (7T7i))lm 



1 _|_ rp i(S-4> 2 ) 



(15) 



X _|_ re «(<5+<M 

where r$ = iV^V^/V^Kdl, r = r ? |r amp |, and 5 ee arg(r amp ) is the relative strong phase. 
The weak phase 2 is defined as the angle between V* b V ud and V t * b V t d, 

v; h v td 



2 = arg . (16) 

\-v: b v ud ) 1 ; 

Unless we know both the magnitude and phase of r amp , it is impossible to extract the 
weak phase 02 from the data a m j X (7T7r). In principle, isospin analysis makes it possible 
to overcome such a pollution without understanding the penguin contribution ||11|| . In 



order to perform the isospin analysis, all modes for B — > tttt have to be measured. 
However, it is difficult to measure the branching ratio of B°(B°) — > 7r°7r°, which has 
background problem as well as a tiny branching ratio of O(10~ 7 ) 0. Therefore, in 
practice, it is hard to perform the isospin analysis. 

In the B — ► K + n~ decay mode, £p = V* b V us is smaller than £p = V t * b V ts by C(A 2 ), and 
the tree contribution does interfere with penguin. For this mode, a m i X (K + 7i~) = 
and the direct CP asymmetry is given by, 

r(g° -» k-*+) - r(g -» K+7T-) 



— 2r sin S sin a;^ 
1 + r 2 — 2r cos 5 cos ' 



(17) 
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where = |^Vt s /V^* b V^ s |. The weak phase u>^ in the sb triangle is defined as ui^ = 
a.Tg(yt s V t * b /V us V* b ) 0. Note that the is related to the weak phase 03 as follows: 

u>£ = * - <j> 3 + - u£, (18) 

where 3 = arg(-V^ V ud /V; b V cd ), w * = a rg(V ud V: s V; d V cs ) = vr + 0(A 4 ). 

We stress that, unless we know the ratio between tree and penguin contribution with the 
relative strong phase, we can predict neither the a d i r nor a m i X . 

Perturbative QCD(PQCD) approach has been developed to theoretically understand 
semi-leptonic and two-body hadronic B meson decays jl2y^,|3| Jl3|Jl4jl . This approach enables 
us to calculate both the phase and magnitude of tree and penguin amplitudes. In this paper, 
applying PQCD approach to B — > nn and B — > Kir decay modes[], we predict the ratios 
r amp between the tree and penguin amplitudes, and give CP asymmetries without relying on 
the isospin analysis. 



II. NUMERICAL RESULTS 

Applying PQCD approach, we take all twist-3 contributions into account, and use the 
wave functions for light mesons, which were decided from light-cone QCD sum rule, 

= ~1^W^{M^ X ) + m ^i x ) + m 07r {fa fa - l)<f%(x)}, (19) 
= ^=f75{/V</4(a:) + m m ct) P K (x) + m m {jb K fa - \)<&{x)}, (20) 

V ^ V c 

where N c = 3 is color's degree of freedom, v£, K s, n^iK) are normalized to dimensionless unit 
vectors, and v^k) ^ P-k(k), n n(K) -L ?V(ir)> v^ac) 'n^ac) — 1- x is momentum fraction of light 
quark's momentum in the meson to parent meson's one. m 0n (K) are defined by the quark 
condensate, 

Ml 2(0|gg|0) Ml (0jgg + S*j0) , . 

m 0n — : — To" > m 0K — j — 72 > \ Zl ) 



where they are given as [IE], 



mon = 1.6 ± 0.2 GeV, m 0K = 1.6 ± 0.2 GeV, (22) 



x In Ref. HH, we used the wave functions for light meson which were obtained phenomenologically. 
In this paper, we consider all twist-3 contributions, and use the wave functions which were decided 
from light-cone QCD sum rule. The updated results hardly change the previous ones [15]. 
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A,P,T 



without SU(3) flavour violation. Lorentz scalar wave functions for light mesons, <r n (K) 
expanded by Gegenbauer polynomials, 



zQx[ 



U ~ x){l - a\ ■ \{\ - 5f) + < • H(i _ 14£ 2 + 21?)}, 

fK Bx(l - x){l - a* -3£-af ■ ^(1 - 5£ 2 )}, 



2y/2N ( 

/K"(7r) J „-K"M 



2 v / 2iV; 



I - • -(1 - 3£ 2 ) + a£ w • -(3 - 30£ 2 + 35£ 4 ) 

/x(tt) 



1 



where £ = 2x — 1. The coefficients a5^r are gi ven as function of m 0K ^ n j 



2x){l±4 W -3(-3±5£ 2 )}, 



)) "2 



are 

(23) 
(24) 

, (25) 
(26) 

and some 



input parameters, 773, ui 3 in Ref. fTEfl, where af , , aj, 773, and 1^3 are calculated from 



QCD sum rule within 30% accuracy. B meson's wave function is parameterized by two 
Lorentz scalar wave functions, 



$ B (x,6) 



1 



V2N C 



\f B + M B ) 75 [0 B (x, b) + ft-)Mx, b)}, 



(27) 



where n + = (1,0, Or), Ti- = (0, 1,0t), and b is the relative separation between 6-quark and 



light quark in the B meson. According to Ref. |13[, the contribution from 0# is found to be 
negligible, and we adopt as the wave function <pB at rest, 



b) = N B x 2 (l — x) 2 exp 



m b x 2 1, ,, 2 
^-2 ( ^ } 



(28) 



where the normalization constant is fixed by the decay constant f B , ub parameterizes 
the extent of B meson, and ujb is order of the mass difference between B meson and b- 
quark, u>b ~ C(A-), A = M B — m b . Note that the decay rate is very sensitive to uj B - Now, 



the branching ratios |Il3||T7| , |T8f will be used to restrict the parameter space. The combined 
branching ratios are 



BR(B° - 


-> 7T + 7T ) 


= (0.44 ±0.09) x 10" 5 , 


BR(B° - 


♦ K+TT-) 


= (1.73 ±0.15) x 10~ 5 , 


BR(B + - 


- K + n°) 


= (1.25 ±0.17) x 10" 5 , 


BR(B + - 


- K\ + ) 


= (1.75 ±0.27) x 10" 5 , 


BR(B° - 




= (1.14 ±0.29) x 10~ 5 , 



(29) 

which are shown in Figure |I[ The errors of the measured branching ratios are still large, and 
it is too early to put them together. Therefore, we obtain the range of allowed parameters so 
that the calculated branching ratios are consistent with the data, eq.(B^), within 2a, where 
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we scan the parameters in the wave function within 0.35 < ujb < 0.55 GeV, = 1.6 ±0.2 
GeV, rriQK = 1-6 ± 0.2, and af, af^\ aj, 773, 003 given in Ref. [I6j within 30% ranges. 

The calculated r amp for each final state are shown at Figs.0, and summarized in Table 
|, where we do not take higher order corrections into account. They show the error for 
I ''amp 1 5 denoted as A|r amp |/|r omp |, is 10 ~ 26%. In principle, a calculation to 0(a 2 s ) should 
be performed. Higher twist contributions should be studied. These will be done in the 
future. For now, we assume that these corrections will not exceed 30% in the amplitude. 
The error A|r amp |/|r amp | is 30 ~ 40% with such higher order corrections. Because the 
factorizable tree amplitude has color suppression in B° — > 7r°7r°, K°tt° decay modes, and 
the tree amplitude comes from only annihilation diagram in B + — > K°tt + decay mode, their 
decay modes have r amp ~ 0(1). Considering KM factor, ~ 2.4 for B — > im and ~ 50 
for B — > Kit, Table | shows that the penguin contribution is dominant in B — > K°ir + , K°7r° 
decay modes and comparable with that of tree in B —>■ 7t tt , K + tt~ , and K + n° decay modes. 
In PQCD, the decay amplitudes are given by integrating the convolution of wave functions, 
hard amplitude, and Wilson's coefficient, over loop momenta. Therefore, the dependence on 
the wave functions does not cancel completely in the ratio r amp . This leads to blobs shown in 
Figs.@. Uncertainty in KM factors, decay constants, and parameters in the wave functions 
exerts an influence on the branching ratios. We elect the central values for the KM factors 
i~9| , |2"0"H and B meson's decay constant /# = 200 ± 30 MeV in Figs.|[ Figs.|2| show the 



value of r amp for each set of parameters which predict the branching ratios to within 2a of 

We prefer not to mix the error from the parameter space and the error from the higher 
order corrections. So, for now, we leave aside the error from the higher order corrections. 
From the allowed range of r amp , we compute the CP asymmetry using eqs.(|2]),(|^). The results 
are shown in Figs|| [|. They are summarized in Table |I[ Now the higher order corrections 
are added to the result shown in Table |T| in quadrature. It results in 40 ~ 47% uncertainty 
for ddir- The measured CP asymmetries have been presented in Ref. p2| , |23"| , where Belle 
Collaboration gives —0.25 < A^K^n'* 1 ) (= adi r (K + 7i~)) < 0.37 at 90% confidence level, and 
BABAR Collaboration gives A K7T (= a dir (K + ir-)) = -0.07 ± 0.08 ± 0.02. They are not 
inconsistent with our predictions. 

Although the relative error of \r amp \ is 10 ~ 26% without higher order corrections, the 
mixing induced CP asymmetry for B — > 7r + 7r~, a m i X (TT + 7r~), has small uncertainty, 6%. This 
is because the magnitude of r amp is small, \r amp \ ~ 0.1. Since a mix has zeroth order term in 
expansion of |r omp |, the error Aa mix /a m i X is proportional to r: 

— ^-(7r + 7r~) = — r sin 02(sin 5 — cos 5 — cot 20 2 sin 5) ■ amp ■ + 0(\r amp \ 2 ). (30) 

fornix I ^amp | 

The factor r = r^|r amp | ~ 0.3 leads to small uncertainty for a m i X . The error from the higher 
order corrections is also suppressed by the factor r, and the total error for a m i X is estimated 
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at 11%. In contradistinction to a m i X , the direct CP asymmetry has large uncertainty. Since 
Acidir/cLdir is not suppressed by r, 

^(tt+O = + <D(\r amp \), (31) 

Qdir | Tamp \ 

the direct CP asymmetry is more sensitive to r amp compared to the mixing induced CP 
asymmetry. Similar errors are assigned to B — > K + 7t~ and K + n°. 

The mixing induced CP asymmetry for B — > Kstt , a mix (Ks^ ), has little uncertainty. 
Here the factorizable tree amplitude is subjected to color suppression, and the penguin 
amplitude without KM factor is comparable to that of the tree, r amp ~ 1. However, KM 
factor, = I^Vtg/V^V^al ~ 50, is very large and it makes the penguin contribution 
dominant. Since r^|r am? ,| 3> 1, the first term in the numerator and denominator can be 
neglected in eq.(|), and A|r amp | is canceled out as shown in section 1. This is the reason 
why a m i X (Ks7T°) has little uncertainty. r^|r am? ,| ~ 50 is classified into the case 3 in section 
1, and a mix (K S 7i°) is given by, 

a mix {K s ^) = sin 2(0! - lu s J), (32) 
where its uncertainty from r amp is 0.03%. 



III. CONCLUSION 

In order to extract the weak phase from the experimental data on CP asymmetries, we 
have to know the magnitude and strong phase of the penguin contribution. We have applied 
PQCD approach to B — > ixix and B — > Kit decay modes, which allows us to compute both 
the strong phase and magnitude of these amplitudes. We have calculated the ratios r amp 
between the penguin and tree amplitude which are crucial in predicting CP asymmetries. 
For each decay mode, uncertainty of \r amp \ is 30 ~ 40%. The direct CP asymmetry which 
is proportional to r amp suffers from uncertainty in |r amp | directly, and the error is estimated 
at 40% ~ 47%. However, the mixing induced CP asymmetry for B — ► 7r + 7r~, which is not 
proportional to r amp , has small uncertainty, and the asymmetry has been predicted with 
11% uncertainty. While it is difficult to measure the time dependence of a^ s7r o at current B 
factories, we predicted r amp for B — > K°7i° decay mode. This leads to the mixing induced 
CP asymmetry which is related to the weak phase <fi\ — iof t . 

There are other attempts to study CP violation in B — » mi, Kix decays [p4 |. In the 
analysis of Ref. [24], SU(3)i? symmetry is parameterized and r amp is obtained from the 



analysis of Ref. ||25||. Our value for r amp which corresponds to the penguin parameter de %e 



(0.24 ~ 0.36) e^ 1010 - 1300 ) in the notation of Ref. @ seems to yield a better agreement with 
experiments. 



S 



Improvement in experimental data is expected in the near future. The uncertainty in 
r amp will be also expected to improve. The predicted CP asymmetries will make it possible 
to extract the weak phase from their experimental data. 
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FIG. 1. The branching ratios for B — > tt + tt and B — > measured at CLEO, Belle, and 
BABAR. The combined data are also shown here. 
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FIG. 2. The ratios 

Tamp between the penguin and tree amplitude for B — > irir, Kir decay modes. 



' amp 

We do not show here r amp for 



7r + 7r° because of \r amp (Tr ' ;~ 



" 1 0. For r amp (7r + 7r ), the value 



obtained here corresponds to de id = (0.24 ~ 0.36) e *( 101 °'" 1300 ) in the notation of Ref. 
should be compared with the value, de %e = 0.09[0.18] e* 193 ! 187 ] ^ given in Ref. 25]. 



This 



13 



1 



T 

cq 



t3 

C3 



0.5 



-0.5 




30 




60 90 120 150 

2 (degree) 



180 



+ 

T 

cq 



0.5 



-0.5 



30 60 90 120 

03 (degree) 



+ 

T 

cq 



0.5 



-0.5 



150 180 



30 60 90 120 150 

3 (degree) 



180 



FIG. 3. Direct CP asymmetry for B — ► tttt, Kit decay modes. The central value of KM factors 
[ 19 , 20 1 gives the darker shaded regions, and the lighter shaded regions include the error of KM 
factors. Orfj r .(7r + 7r°) is almost zero for any 02- a dir{K°7r + ) is almost zero for any 03. a^j r (i^°7r°) 
becomes maximum at 03 = 90°, adi r {K°TT°) = —0.04. 
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FIG. 4. Mixing induced CP asymmetry for B — > ir + Tr~. The difference from the dotted 
line(sin 202) shows sizeable penguin pollution. 
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TABLES 



final state \r amp \ arg(r amp ) 

vr+TT- 0.089 ~ 0.134 108° ~ 135° 

7r vr 1.11 ~ 1.43 -144° 130° 

vr+vr 0.0113 ~ 0.0146 180° 

K + ir- 0.087 ~ 0.106 137° ~ 155° 

K + tt° 0.095 ~ 0.117 143° ~ 159° 

K°ir+ 1.91 ~ 2.85 112° ~ 128° 

K°tt° 0.95 ~ 1.60 92° ~ 122° 



TABLE I. 

r amp f° r each final state calculated in PQCD. Uncertainty conies from parameters 
in the wave functions. 



final state 


4>i 




a dir 




7T + 7T~ 


80° 


0.366 


~ 0.550 (0.298 ~ 0.636) 


20% (36%) 


7T+7T 


any 




~ 




K+TT- 


70° 


-0.257 ~ 


-0.180 (-0.261 ~ -0.153) 


18% (26%) 


K + TT° 


70° 


-0.205 ~ 


-0.130 (-0.211 ~ -0.108) 


22% (32%) 


K TT + 


any 




~ 




K°7T° 


90° 


-0.0393 ~ 


-0.0237 (-0.0428 0.0213) 


25% (33%) 




02 




Q-mix 


^&mix/ fornix 


TT + TT~ 


40° 


0.904 


~ 0.966 (0.864 ~ 0.979) 


3.3% (6.2%) 



TABLE II. For illustration purposes, we show direct CP asymmetries and mixing induced CP 
asymmetry for B — > tttt, Kit decay modes. We have chosen fa (fa = fa for B — ► tttt and fa = fa 
for B — > Kit) which give the maximum asymmetries. The parenthesized values include errors of 
KM factors, while the unparenthesized values are obtained from the central values of KM factors 
[ 19 , 20 . The higher order corrections are not considered in the above errors. 
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